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Abstract. We develop a quantum mechanical formalism to treat the strong coupling 
between an electromagnetic mode and a vibrational excitation of an ensemble of organic 
molecules. By employing a Bloch-Redfield-Wangsness approach, we show that the 
influence of dephasing-type interactions, i.e., elastic collisions with a background bath 
of phonons, critically depends on the nature of the bath modes. In particular, for 
long-range phonons corresponding to a common bath, the dynamics of the “bright state” 
(the collective superposition of molecular vibrations coupling to the cavity mode) is 
effectively decoupled from other system eigenstates. For the case of independent baths 
(or short-range phonons), incoherent energy transfer occurs between the bright state 
and the uncoupled dark states. However, these processes are suppressed when the 
Rabi splitting is larger than the frequency range of the bath modes, as achieved in a 
recent experiment [Shalabney et al ., Nat. Commun. 6, 5981 (2015)]. In both cases, 
the dynamics can thus be described through a single collective oscillator coupled to a 
photonic mode, making this system an ideal candidate to explore cavity optomechanics 
at room temperature. 


1. Introduction 

The field of cavity optomechanics explores the interaction between electromagnetic 
radiation and the quantized mechanical motion of nano- or micro-oscillators [1-3]. 
Recent developments promise a rich variety of applications such as precision mechanical 
measurements and coherent control of quantum states (see [4, 5] for recent reviews). 
When the interaction becomes sufficiently large, coherent energy exchange between the 
optical and mechanical degrees of freedom becomes possible and the system could reach 
the strong coupling regime. Along this direction, a novel approach consists in coupling a 
cavity mode to a molecular bond vibration [6, 7], which is in the ground state already at 
room temperature, without requiring any cooling. In order to achieve strong coupling, 
this has to be done by using an ensemble of molecules. This phenomenon has common 
ingredients with the collective strong coupling observed when electronic excitations of 
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Figure 1 . Sketch of the model system. An ensemble of N molecules interacts 
collectively with a cavity mode of frequency cj c . The cavity is tuned to a molecular 
vibrational mode. The remaining rovibrational modes are introduced by either 
connecting all molecules to a common bath (panel a) or by coupling each molecule with 
an independent bath (panel b). In both cases the harmonic bath is characterized by a 
spectral density J(w). 

organic molecules interact with cavity [8-10] or plasmonic modes [11], which has been 
studied extensively during the last years (see [12, 13] for recent reviews). 

In the strong coupling with an ensemble of vibrational modes, the cavity resonance 
couples to a collective superposition of the molecular vibrations, the so-called bright 
state, forming hybrid states called polaritons. In principle, all other superpositions of 
vibrational excitations (the so-called dark states) remain uncoupled. However, within 
this setting, a thermal bath of low-frequency rovibrational modes, which normally only 
introduces dephasing [14, 15], interacts with the vibrational excitations and may influence 
the system dynamics. It is an open question how the dephasing mechanisms affect the 
strongly coupled dynamics and what role the dark states play. In particular, it is 
unknown whether the bright mode in such a system effectively behaves like a single 
isolated oscillator, which would enable the direct application of protocols developed in 
cavity optomechanics. 

In this paper, we present a fully quantum theory of the phenomenon of collective 
strong coupling of molecular vibrations with a cavity mode. In contrast to classical 
transfer matrix calculations, which can be used to fit the experimental spectra [6, 7] but 
only provide very limited information about the incoherent dynamics in the system, our 
approach naturally incorporates all incoherent processes, in particular those induced by 
dephasing-type interactions. We describe the system using a quantum-mechanical model 
of a single photon mode coupled to an ensemble of harmonic oscillators representing 
molecular vibrations (the C=0 bond stretching mode of polyvinyl acetate at an energy 
of 215 meV in the experiment [6]). In order to incorporate the coupling of the molecular 
vibrations to low-frequency rovibrational modes, we assume that the oscillators are 
connected to either a common or to independent phononic baths (see Fig. 1). We employ 
Bloch-Redkeld-Wangsness (BRW) theory [16, 17] to obtain master equations describing 
the system dynamics under strong coupling. The final master equations only contain a 
few Lindblad terms, with rates determined by the product of i) the bath spectral density 
at the corresponding transition frequency and ii) algebraic prefactors obtained from 
transforming the system into its eigenstate basis. Using the specific properties of the 
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system under study, we can evaluate all these prefactors analytically , and are thus able to 
directly read off the population transfer rates between the system eigenstates. This allows 
detailed insight into the system dynamics and, specifically, the role of the dark states. 
Our results demonstrate that for large enough Rabi splitting, the bright state indeed 
behaves like a single isolated oscillator and the dark states play an almost negligible role. 
For the case of a common bath, corresponding to long-range bath phonons, this is even 
true regardless of the Rabi splitting. 


2. Theory 


2.1. Coherent dynamics 

We model the system as a set of N molecular vibrational modes coupled to a single 
electromagnetic (EM) mode in a microcavity, as depicted in Fig. 1. The mirrors in the 
experiments [6, 7] are actually planar, and the photonic modes form a continuum, with a 
dispersion relation depending on the in-plane momentum k. However, the assumption of 
a single EM mode coupling to many molecules is justified when comparing the density of 
EM modes with the molecule density in the experiment. The 2D density of EM modes 
with wavevector k < k max (where k = |/c|) is N p h = J 0 fcmax 2nkdk = x /(47t). For 

the jth transversal mode in a cavity of length L with background refractive index n, 
the dispersion relation is uj = (^-) 2 , so that the density of states with energy 

uj < uj m ax is (here and in the following, we set h = 1 ) 



_rf_ ( 2 

4t rc 2 ^ max 



(1) 


where ujj = fgj-. We assume that modes with j — 1 and uj < u; max ~ oj\ + Dr participate 
in the dynamics. Taking the parameters from the experiment of Shalabney et al. [6] 
(n ~ 1.41, uj\ = 215 meV, Dr = 20.7 meV, L ~ 2 p m), this leads to an EM mode density 
of IVph 4 • 10 6 cm -2 . In contrast, the reported molecular density is d ~ 8 • 10 21 cm -3 , 
giving a 2D molecular density of Ld ~ 2 • 10 18 cm -2 . There are thus on the order of 
10 12 molecular vibrational modes coupled to each photonic mode, and our assumption is 
clearly justified. 

Within the rotating wave approximation, i.e., neglecting processes that create or 
destroy two excitations, the coherent dynamics of the system is then governed by the 
Hamiltonian: 

N N 

H s = uo c a)a + ^2 + XI 9i( ac l + a ^ c i), (2) 

2=1 2=1 

where a is the annihilation operator for the cavity mode with frequency uo c , and c % is the 
annihilation operator of the optically active vibrational mode of molecule i, characterized 
by its frequency uj 1 and its position ry The cavity-oscillator interaction is given by g t , 
which depends on the cavity electric field strength and the change of the molecular 
dipole moment under displacement from the equilibrium position (see [6] for details). 
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For simplicity, in this work we assume a regular configuration in which all the molecules 
are identical (cn* = cu m , gi = g), as well as zero detuning (u; c = uj m ). We have explicitly 
checked that orientational disorder (/ gf) and inhomogenous broadening (07 ^ uf) do 
not significantly affect the results presented below. 

We focus on the linear response of the system, so that we can restrict the treatment 
to the zero- and single-excitation subspacesfi The N + 1 singly excited eigenstates of H s 
are formed by: i) two polaritons, |±), symmetric and antisymmetric linear combinations 
of the cavity mode a^O), with the collective bright state of the molecular excitation, 
| B) = |±) = ;^(g^|0) ± | B)) and ii) the so-called dark states, N — 1 

combinations of molecular excitations orthogonal to \B), which have eigenfrequencies u> m 
and no mixing with the photonic mode. The eigenfrequencies of the two polariton modes 
are u; m ± gV~N, with Rabi splitting fl R = 2 gVN . Collective strong coupling emerges 
when fl R is larger than the losses of the system. We can distinguish three types of loss 
mechanisms: cavity losses (rate k), non-radiative internal losses within the molecules 
(rate 7 nr ), and dephasing-type interactions. Spontaneous radiative decay is very slow 
(on the scale of milliseconds) due to the low transition frequencies and can be safely 
neglected. 

2.2. Incoherent dynamics 

We now turn to the description of the incoherent dynamics induced by the different loss 
mechanisms. Whereas cavity losses and nonradiative internal molecular decay can be 
seen as pure decay channels and included as constant Lindblad terms when analyzing 
the system dynamics, dephasing-type interactions must be treated in a more detailed 
fashion, as we detail below. In the case analyzed in this work, these interactions are due 
to elastic scattering of low-frequency rovibrational bath modes with the main vibrational 
modes involved in strong coupling, described by the interaction Hamiltonian 

N 

H<t> = ^2 c \ c i ^2 ^ ik (pik + b\ k ) . (3) 

i =1 k 

The spatial extension or localization of the bath modes determines the character of the 
coupling. In the following, we focus on two limiting scenarios. In the first scenario, all 
the vibrational modes are coupled to the same common bath (see Fig. la), characterized 
by delocalized phonons bik = bk with bath Hamiltonian Hff ia = oj kb\bk- In the 

second scenario, each molecular vibrational mode is coupled to an independent bath 
characterized by on-site phonons b^, with bath Hamiltonian H]f d = 'fQ lk 
both scenarios, we assume that all molecules are identical (A^ = A k and 17 ^ = ay). 

The frequently used approach of treating dephasing through a frequency-independent 
Lindblad superoperator acting only on the vibrational modes is not valid in our case as, 

f We note that since the system Hamiltonian is quadratic in the bosonic modes {a, c,}, it can in principle 
be diagonalized without restriction to any excitation subspace. This leaves most of the derivation 
presented in the following unchanged, but introduces additional prefactors depending on the excitation 
numbers in the incoherent rates derived below. For simplicity, we thus allow at most one excitation. 
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within the strong coupling regime, the molecule-cavity coupling (frequency Qr) is much 
faster than the correlation time of the phononic environment. Instead, the influence of 
the background modes has to be taken into account in the dressed basis obtained after 
diagonalizing the strong-coupling interaction. The bath is completely characterized by 
the spectral density J*( uj) = A ? k 5(oj — u> ik ) (with Ji(co) = J{yj) in our case). In the 

following, we assume an Ohmic environment with a quadratic cutoff at frequency oj cut 

J(co) = rjcoe~ ( ' UJ ^ UJcut ' >2 , (4) 

where rj is a dimensionless constant that determines the system-bath coupling strength. 

If the system-bath coupling is sufficiently weak, Bloch-Redfield-Wangsness (BRW) 
theory [16, 17] can be safely applied to derive a master equation for the system 
dynamics [18, 19]. This approach relies on the first and second Born approximation, which 
calculate the effect of the system-bath coupling up to second order perturbation theory 
and assume that the bath state remains unmodified, i.e., the bath density matrix pi, is 
time-independent (and thermally populated in the following), where Ptotai(^) — pit) 0 
Additionally, since the decay of the bath correlations ~ 1 /o; cut occurs on a much shorter 
time scale than the dynamics caused by the interaction with the bath, the Markov 
approximation is used, disregarding all memory effects of the system-bath interaction. 
In the interaction picture (denoted by a tilde, O = e ^ Hs+Hb ^ t O e ~ l( ' Hs+Hb ' >t ), ^hc S y S tem 
density operator p then evolves according to 

t 

d t p(t) = - I Tr b [Ht(t),[Ht(lf),m®Pb]W, (5) 

— OO 

where Try denotes the trace over the bath degrees of freedom. 

The bath-dependent part of the system-bath coupling is fully encoded in the 
bath correlation functions between the modes on molecular sites i and j, (f>ij(r) = 
Trb\pl k (T)bjk(0)pb]. It is independent of both i and j for a common bath, 0^° m (r) = 
6(t), while for independent baths, the off-diagonal terms vanish, (/>™ d (r) = The 

autocorrelation function 0 (r) can be expressed through the spectral density [ 20 ], 

OO 

4>(r ) = J J(co) |[n(o;) + l]e~ luJT + n(o;)e* a,r } da;, ( 6 ) 

o 

where n(uj) = (e^ — l ) -1 is the Bose occupation factor and (3 = 1/ksT, with ks the 
Boltzmann constant and T the temperature. 

We now evaluate Eq. (5) within the system eigenbasis. Before we proceed, we 
note that the common approach of using frequency-independent Lindblad terms to 
describe dephasing is equivalent to neglecting the strong coupling in the incoherent 
dynamics, i.e., to use only the uncoupled system Hamiltonian uj c a)a + uj m J0 cjc, in the 
interaction picture on the right-hand side of Eq. (5). When the molecule-cavity coupling 
is comparable to or faster than the decay of bath correlations {BLr > o; cut ), this is an 
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invalid approximation, and it is crucial to include the full system Hamiltonian when 
deriving the master equation to satisfy detailed balance [ 21 ]. 

We thus proceed by expressing the system part of (oc cjc*) in terms of the dressed 
eigenbasis and inserting this expansion in Eq. (5). This leads to 

N °® 

d t p=Y2^2 I 0ij( T ) u ip u iq u jrUj S e i{ “ pq ~ Usr)t+l “ srT [a rs p(t),cr P q]dT + R.C., ( 7 ) 

i,j=l p,q,r,s J Q 

where a a b = |a)(6|, H s \a) = 0J a \a), uj a b = cu a — u>b, and the sums over p , q, r and s include 
all system eigenstates. Furthermore, Ui a = (a|c|| 0 ) give the overlaps between system 
eigenstates and vibrational mode excitations and can be chosen real. Finally, we made 
the substitution t' — t — t. Inserting ffr) from Eq. ( 6 ) leads to integrals of the type 

e ±iAur dT = ^(Au;) ± i p.v.(Acc -1 ), (8) 

where P.V. denotes the Cauchy principal value—we neglect these imaginary parts as 
they only induce small energy shifts ( Lamb shifts ) that can be reabsorbed in the coherent 
dynamics, and arrive to 



N 

d t p=Yl ^2 Sij(Usr)u ip U iq U jr U js {e l (“ pq ~ Usr)t [i cr rs p(t ), <J pq \ + H.c.} , (9) 

i,j =1 p,q,r,s 

with the bath noise-power spectrum 


S(u>) 


TrJ(cu)[n(uj) + 1 ] oj > 0 
7T J (— £n)u( — £n) OJ < 0 


( 10 ) 


For a common bath, = S{uj), while for independent baths, Sf- A (uj) = S ltj S(cu). 

Note that ^fO) = lim S(co) = lim S(co) is well-dehned if lim exists. 

In the resulting master equation, the terms in which u pq differs from uj sr oscillate 
as a function of time t, which can lead to a violation of the positivity of p for long 
times. If the timescale of the bath-induced system dynamics is much slower than 
the coherent dynamics, i.e., r ( p rsc, these terms can be removed by averaging the 
master equation over a time short compared to r^, but long compared to rsc [22, 23]. We 
can estimate Tp, ~ l/q^, with 7 ^ the bare-molecule dephasing rate, for which the total 
bare-molecule width presents an upper bound, 7* < 7 = 3.2 meV [ 6 ]. As a consequence, 
only the secular terms uj pq = uj sr persist. This secular approximation 1 results in a master 
equation where populations and coherences are decoupled. We note that the secular 
approximation aids in the interpretation of the different terms that are obtained, but is 
not required and indeed only used in some of the results shown in the following. 

f The secular approximation is also sometimes called the rotating wave approximation , but should not 
be confused with the rotating wave approximation performed in the system Hamiltonian in Eq. (2). 
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Table 1. Secular terms in the master equation and the labels for the associated “bare” 
rates g rs = 2 S(w sr ). The dark states are labeled with d £ 1,2,JV — 1. 

The secular terms are enumerated in Table 1, which lists the states { pq , rs} connected 
by transitions with frequency u sr , as well as the labels we assign to the “bare” rates 
7 rs = 2S(u sr ), which do not contain the algebraic prefactors. The definition is chosen so 
that for terms where pq = sr, we obtain 

S(pi S r)\\(JrsPiff i &si\ T H.C.} 7rs-^<r rs [p] j (H) 

where Cx[p\ = XpX1 | {X^X, p} is a standard Lindblad superoperator. Positive 

frequencies uj sr > 0 correspond to phonon emission where the system transitions from a 
higher- to a lower-energy state, while negative frequencies cu sr < 0 correspond to phonon 
absorption. We obtain secular terms connecting the two polaritons (% = 25(Hr), 
7 a = 25(— Hr)§), terms connecting the polaritons with the dark modes (T e = 25(Hr/2), 
T a = 25(— Hr/2)), and terms connecting states with the same energy ( 7 ^ = 25(0), equal 
to the bare-molecule dephasing rate). The latter give pure dephasing for the polaritons 
|+), |—), but produce coupling between populations and coherences for all dark states. 

The final master equations are obtained by using the properties of the 
basis transformation matrix u ia to evaluate the algebraic prefactors JT Ui p Ui q Ui r Ui s 
(independent baths) and Yhi j u ip u iq u jr u js (common bath) in Eq. (9). Specifically, we 
use that i) the polariton-vibrational mode overlaps are given by = ±^=, ii) that 
dark states are orthogonal to each other QTk u ui u uV = ( fu'), and iii) that dark states are 
orthogonal to the polaritons (Xr^W = 0). Note that i and j in the sums are molecule 
indices (not including the cavity mode), so that Yli u ip u iq = &p,q is generally not true. 

This procedure gives the final secularized master equation for the density operator, 
given in the Schrodinger picture below. For a common bath with 4>ij(r) = 5 (r), we find 
that many terms vanish because of the orthogonality relations, giving 

d t p = -i[H s , p] + C a+ _ [p] + C a -+ [p] + ^ £< t pp [p] + 7 '<t>£v[p\, ( 12 ) 

p=+- 


The same 7 a ,7e are obtained under strong classical driving of a single two-level system [24]. 
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Figure 2. Illustration of the different decay and dephasing mechanisms that emerge 
from our modeling for the two dephasing scenarios: common bath (panel a) or 
independent baths (panel b). The arrows indicate decay processes while wavy lines 
represent elastic events that produce dephasing. The fuzzy halo around the dark states 
indicates excitation transfer and dephasing acting within this manifold. 


where V = add is the projector into the dark-state subspace. The Lindblad terms 
Cx[p] correspond to incoherent excitation transfer between system eigenstates and are 
depicted schematically in Fig. 2a. A phonon of frequency Qr may be emitted transferring 
an excitation from the upper polariton |+) to the lower polariton |—), with rate 7 e /4. 
Phonon absorption occurs analogously, with characteristic rate 7 a /4. Furthermore, the 
polaritons undergo pure dephasing with rate 7^/4. The last term corresponds to bare- 
molecule dephasing for a common bath, but projected into the degenerate dark-state 
subspace (using 'E ) Yli c \ c iD = T). Remarkably, for the case of a common bath, i.e., 
long-range bath phonons, the dark states are completely decoupled from the polaritons 
and the bright state behaves identically to a single oscillator interacting with the cavity 
field. 

Turning to the case of independent baths, Aj( T ) — Aj ( M r )> we instead find 


d,P — —i[H s , p] + ~^C a+ _ [p] + -^C (T - + [p] 

(13a) 

+ M + £<r+d\p\) + 

d d 

->]) (13b) 

p _ ^ 

+ 4^ (\- a d-P’ a d+] - [<T+dP , <r-d] + H.C.) 

d 

(13c) 

P _ 

+ a d-} - \ a -dP1 v+d] + H.C.) 

d 

(13d) 

+ 4N ^2 + 70 ^2 t-"Dc\aM > 

P=+~ i 

(13c) 

which now includes excitation transfer between the polaritons and dark states (r a , 
r e ), driven by phonons of frequency Qr/2. While the individual terms are strongly 
suppressed by the prefactor 1/N (with N ~ 10 12 in the experiments), excitation transfer 
into the dark states still occurs efficiently due to the sum over N — 1 dark states d. For 
N — y oo, the total population transfer from the polaritons to the dark states thus occurs 
with rates r a /2 and r e /2, as shown in Fig. 2b. On the other hand, pure dephasing 


of the polaritons and direct transitions between them through absorption or emission 
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of phonons of frequency Qr play a negligible role within this dephasing scenario, as 
their rates scale as 1/N. Additional terms couple between different polariton-dark state 
coherences (Eqs. (13c), (13d)), without affecting the populations. Finally, the second 
term in Eq. (13e) again corresponds to a bare-molecule Lindblad dephasing term that 
has been restricted to act only within the degenerate dark-state subspace. For N —» oo, 
Eq. (13) simplifies to 

dtP = —i\H s , p} + -y £v-[P\ + yA?+[p] + 7 <t> ^ ^v c \av\.P\ > (14) 

i 

where jCr± = d C ad± is an averaged Lindblad superoperator inducing equal 

population transfer from a polariton to all dark states. For large N, the dark states thus 
act like a sink, and any population transferred to them is trapped and does not further 
participate in the polariton dynamics. 

As is clear from the expressions for the different decay rates, Qr is the parameter 
that controls the interaction between the dressed excitation and the phonon bath. Hence 
the importance of these decoherence mechanisms is dictated by the spectral density 
evaluated at f Ir or VLr/2. This is in clear contrast to what we would obtain by treating 
the interaction H 0 in the uncoupled basis. In that case, the Markov approximation would 
have resulted in “standard” Lindblad terms 7 t (p) (common bath) or 7 ^ JA C t ip) 
(independent baths). In both cases, the dephasing interaction would be totally controlled 
by just the zero-frequency limit of the spectral density, resulting in an overestimation of 
the rates 7 ^ and r a e . This fact emphasizes the key importance of deriving Lindblad 
terms in the strongly coupled basis when considering non-flat reservoirs [ 21 ]. 

Apart from the decay and dephasing mechanisms induced by the dephasing-type 
interactions, which conserve the number of excitations, the excited states may decay 
through nonradiative molecular decay and cavity losses into the ground state |0). These 
decay process, which we have neglected in the theory up to now, could be described 
within the same framework by dissipative coupling with thermal baths (e.g., the photonic 
modes outside the cavity). However, as the energy shifts induced by the strong coupling 
are small compared to the transition frequency (Qr <C 0J m ) , we can include them through 
Lindblad terms for the bare cavity and molecular vibrational modes. After removing 
nonsecular terms in the eigenstate basis, this gives new Lindblad terms with decay rates 
°f (7nr/2 + k/2) for the two polaritons and 7 nr for the dark states. 

3. Results 

In the following, we apply our theoretical framework to the experimental results of 
Shalabney et al. [ 6 ]. For the vibrational mode of the bare molecules, they report a 
linewidth of 7 = 3.2 meV, with negligible inhomogeneous broadening. This linewidth 
has contributions from nonradiative decay and dephasing, 7 = 7 nr + 7 ^, which can not 
be distinguished in the absorption spectrum. Although direct information about the 
weights of non-radiative and dephasing channels is thus not available, dephasing is in 
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uj [meV] uj [meV] 

Figure 3. Absorption spectra for the two possible dephasing scenarios as a function 
of /. In both panels the two sets of curves are offset for a better visualization. In 
panel (a), the parameters used for the calculations are those reported in Ref. [6] (see 
text), whereas in panel (b), the Rabi splitting and the cavity decay rate are reduced to 
f Ir = 6.5 meV and k = 1 meV, respectively. 


general expected to provide a significant contribution for vibrational transitions [14, 15]. 
Therefore, in our calculations we will use a factor / [ 7 ^ = /q, q nr = (1 — /)q] to 
measure the relative importance of the two channels. In this way, the factor 7 in 
Eq. (4), which quantifies the strength of the system-bath coupling, is simply given by 
rj = /y/(27 t/ 1 , /3To) where To = 300 K. The cut-off frequency for the thermal bath of 
low-frequency rovibrational excitations is chosen as oy ut = 6 meV, corresponding to the 
range of low-frequency phonon modes in the system [ 6 ]. Finally, we use a cavity loss 
rate of k = 17 meV as estimated in [ 6 ] through fitting of the transmission spectrum in 
the strong coupling regime. 

We calculate the absorption spectra by introducing a weak driving term, Hjft) ~ 
ae~ luJt + H.c., which coherently pumps the cavity mode. The density matrix in the 
steady state, p ss , can be calculated in the frame rotating with the driving frequency 
u. In this frame, H ( j is time-independent, but the system frequencies are shifted and 
the density matrix p ss depends on uj. The absorption spectrum is then obtained as 
A(uj) oc Tr[p ss (a;)a]. For the numerical implementation, we employ the open-source 
QuTiP package [25]. 

Fig. 3a depicts the theoretical absorption spectra for the parameters reproducing 
the experimental situation, with Rabi splitting fl R = 20.7 meV, for the two possible 
dephasing scenarios analyzed in this work (common or independent baths). Different 
values of / are tested: / = 0, / = 0.5 and f — 1. Both approaches coincide in the limit 
/ = 0 (no dephasing), but their behaviour differs for non-zero /, as can be inferred from 
the different decay rates as rendered in Fig. 2. The Rabi splitting is much larger than the 
range of low-frequency vibrations w cut ) and, hence, all terms connecting levels 

with different energies ( 7 a , q e , F a , T e ) are essentially zero. The widths of the polaritons in 
the absorption spectrum (for N —y 00) are simply given by q± = K+ ^" r = | + (1 — |) | 

for a common bath, and q± = K+ 2 ynr = | + (1 — /)| for independent baths. 

This also implies that for both bath scenarios and within the experimental conditions 
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Figure 4. Population dynamics starting from the initially excited upper polariton |+) 
for the two possible dephasing scenarios: (a) Common bath and (b) independent baths. 
The parameters are the same as in Fig. 3(b), with / = 0.5 and N = 100. Each panel 
shows the upper polariton (+) and lower polariton (—) populations, as well as the total 
population in all dark states (DS), given by Tr (Dp). 


reported in [6, 7], the ensemble of vibrational modes behaves effectively as just a single 
collective molecular oscillator, the bright state, coupled to the cavity mode. Consequently, 
the dark modes are effectively decoupled from the system dynamics under external driving. 
We note here that the fit used to extract the cavity width k from the experimental 
data in [6] is performed under strong coupling (using a transfer matrix method). Thus, 
the change of linewidths predicted by our model would already be present in the 
observed spectra, and is consequently absorbed in the extracted cavity linewidth. This 
unfortunately precludes a direct test of our model based on the experimental absorption 
data, for which the bare cavity linewidth would need to be available. 

Furthermore, it is interesting to note that in the case of individual baths, the 
dephasing contributions to the polariton modes are completely suppressed, analogous 
to the well-known suppression of inhomogeneous broadening under strong coupling [26]. 
This can be understood by the fact that dephasing and inhomogeneous broadening 
are closely related, corresponding respectively to temporal fluctuations or to a static 
distribution of the transition energies. 

The dark modes could play a bigger role in a situation with smaller Rabi splitting, 
for which, in order to still achieve strong coupling, the cavity losses also need to be 
reduced as compared to the experiments (e.g., by using thicker mirrors). The resulting 
absorption spectra are shown in Fig. 3b, for the same parameters as in Fig. 3a but now 
with k = 1 meV and Qn = 6.5 meV. In both bath scenarios, a slight asymmetry between 
upper and lower polariton is now noticeable, as the emission of phonons from the upper 
polariton (oc n[uf) +1) is more likely than the absorption of phonons in the lower polariton 
(oc n(u)). However, as the involved transition frequencies are smaller than the thermal 
energy (UbT = 25.9 meV for T = 300 K), the thermal occupation n{uS) is significant and 
the rates for phonon emission and absorption are comparable. Furthermore, although 
the phenomenology in the observed absorption spectra is quite similar in both dephasing 
scenarios, with a reduction of the polariton linewidths for increasing dephasing /, the 
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underlying physics are now quite distinct. This is demonstrated clearly when inspecting 
the population dynamics, shown in Fig. 4 for the case of the upper polariton being 
initially excited. For the common bath, shown in Fig. 4(a), the dark states are completely 
decoupled from the dynamics and population is only transferred between the polaritons. 
In contrast, in the case of independent baths, shown in Fig. 4(b), the population primarily 
is transferred from the upper polariton to the dark states, which then decay through 
nonradiative losses. 

4. Conclusions 

In summary, by using a fully quantum framework we have studied in detail the 
phenomenon of collective strong coupling when an ensemble of molecular vibrational 
modes interacts with a cavity electromagnetic mode, as realized experimentally in two 
recent papers [6, 7]. We have demonstrated that dephasing-type interactions with a 
thermal bath of background modes in such a system have to be treated beyond the 
usual Lindblad approximation in order to represent the effects of the spectral density of 
bath modes correctly. We have investigated two “extreme” scenarios for the bath, with 
either a common bath for all molecules, or independent baths for each molecule. For 
the experimentally relevant parameters, we find that the dark modes are almost totally 
decoupled from the polaritons in both scenarios, and the bright state behaves almost 
like a single isolated oscillator. Our findings thus suggest that this type of system is an 
ideal and simple platform to explore the exciting possibilities of cavity optomechanics at 
room temperature. 
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